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THERMAL  AND  CREEP  EFFECTS  IN  WORK-HARDENING  ELASTIC-PLASTIC 

C<  AT  TAA  1 

By  Alexander  J.  Wang  and  William  Prager  (Brown  University) 

Extremum  principles  governing  the  isothermal  deforma- 
tion of  a work-hardening  elastic-plastic  solid  have  been  given 
by  Hodge  and  Prager^  and  Hill^.  In  the  present  note  it  is  shown 
how  these  principles  can  be  extended  to  include  thermal  and  creep 
effects. 

Using  rectangular  Cartesian  coordinates  x^  (i  = 1,2,3), 
denote  the  infinitesimal  displacement  from  the  standard  state  by 
ui?  the  infinitesimal  strain  by  e ^ ^ , tho  stress  by  d ^ ^ , the  mean 
normal  stress  by  o,  and  the  stress  deviation  by  s-y.  The  mean 
normal  stress  is  defined  as 

d = 3 dii  > 

where  the  usual  summation  convention  regarding  repeated  subscripts 
is  used;  the  stress  deviation  is  defined  as 

si j = dij  “ d 6ij  » 


*The  results  presented  in  this  paper  were  obtained  in  the  course 
of  research  conducted  under  Contract  N7onr-35801  between  the 
Office  of  Naval  Research  and  Brown  University. 

2P.  Hodge  and  W.  Prager,  A Variational  Principle  for  Plastic 
Materials  with  Strain  Hardening,  jrMath.Phys.  ~?7.  l-loCT^S). 

^R.  Hill,  The  Mathematical  Theory  of  Plasticity.  Clarendon  Press, 
Oxford,  1950,  pp. 63-66. 
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6.jj  being  the  Kronecker  delta.  The  mean  normal  stress  is  an 
invariant  cf  the  stress  tensor.  Further  invariants  useful  in 
the  theory  of  isotropic  plastic  solids  are 

J2  = \ sij  sji  » J3  = | sij  sjk  ski  * (3) 

To  specify  the  mechanical  behavior  of  the  isotropic 
solid,  consider  an  element  that,  at  the  generic  instant  t,  has 
the  temperature  9 and  the  strain  and  is  under  the  stress  c^y 
In  the  interval  between  the  instants  t and  t + dt,  let  the 
temperature  be  changed  by  d©  and  the  stress  by  dd^  = ds^ 

+ dd&^j . 

The  corresponding  change  of  strain,  de.^,  will  then  be 
assumed  to  consist  of  the  following  components: 

1)  the  elastic  component 

de^  = a(0)dSlj  + P(9)ddbij  , (4-) 

where  a(Q)  is  one  half  of  the  reciprocal  of  the  terperature- 
dependent  shear  modulus,  and  p ( 0)  is  one  third  of  the  recip- 
rocal of  the  temperature-dependent  bulk  modulus; 

2)  the  thermal  component 

de^j  = [a‘(0)  s±i  + P'  <9)d&1;j  ] dO  + Y(Wi>}  <*©  (5) 

where  the  prime  denotes  differentiation  with  respect  to  the 
temperature  and  y(9)  is  the  coefficient  of  linear  thermal  ex- 
pansion at  zero  stress; 

3)  the  creep  component 
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where  9 = cpCO,^,^);  and 
4)  the  plastic  component 


3 

(6) 


0 if  ^(e,J2,J3)  < hj 

T(©,J2,J3)|^ — (dij)  + 1 d\Ji|  ) if  4)(G,J2,J3)  = h, 


(7) 


where  the  scalar  h describes  the  state  of  hardening  of  the 
considered  element  at  the  time  t. 

Equations  (4)  and  (5)  result  from  a generalized  form  of 
Hooke's  law  in  which  the  elastic  constants  are  functions  ol  the 
temperature. 

In  Eq.(6)  the  derivative  a<p/Q<5  must  be  evaluated  for 
the  state  of  stress  existing  at  the  time  t.  The  expression 


4>  9cp/0cj  then  represents  the 
ing  to  this  state  of  stress. 


rate  of  secondary  creep  correspond- 
This  expression  is  sufficiently 


. general  to  comprise  all  isotropic  laws  of  secondary  creep  that 


have  so  far  been  proposed  in  the  literature  on  creep. 

In  (7),  the  equation  \K©,J2,J3)  = h represents  the 
temperature-dependent  yield  limit  for  the  state  of  hardening 
achieved  at  the  instant  t.  The  first  line  of  (7)  therefore  states 


that  there  will  be  no  change  in  plastic  strain  where  the  state  of 
stress  is  below  the  yield  limit.  The  differential  d^  in  the 
second  line  of  (7)  must  bo  evaluated  from  the  given  temperature 
and  stress  at  the  instant  t and  the  changes  in  temperature  and 
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stress  during  the  considered  time  interval.  On  account  of  the 

absolute  value  in  the  second  line  of  (7),  there  will  bo  no  change 

in  plastic  strain  even  if  the  state  of  stress  is  at  the  yield 

limit  provided  that  d4>  < 0. 

The  extremum  principles  that  are  to  bo  established 

concern  the  following  boundary  value  problem.  Consider  a mass 

of  work-hardening  plastic  material  that  has  been  deformed  and, 

at  the  time  t,  occupies  a region  V bounded  by  the  surface  S. 

Suppose  that  the  temperature  0,  the  stress  d . . , and  the  state 

~ J 

of  hardening  h are  known  throughout  V.  If  the  unit  vector  along 
the  exterior  normal  of  S is  denoted  by  , the  surface  traction 
in  the  considered  state  is  = d^n^.  Prescribe  now  infinites- 
imal changes  d©  of  the  temperature  throughout  V,  infinitesimal 
changes  dT.^  of  the  surface  traction  on  the  portion  ST  of  the 
surface  and  infinitesimal  displacements  dU^  on  the  remainder 
Sy  of  the  surface.  What  are  the  corresponding  changes  of  stress 
doij  and  the  corresponding  displacement  du.^  throughout  V? 

The  infinitesimal  displacement  du.j_  causes  the  strain  to 

change  by 

delJ  “ $ CjfrWuJ)  + Ux7^aui^  ’ (8) 


The  change  of  stress  must  satisfy  the  equation  of  equilibrium 
which,  to  within  higher  order  terms,  can  be  written  as 

" »•  « 
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Finally,  the  following  boundary  conditions  must  bo  satisfied 


ddij  n^  = dT.j_  on  ST  , 


du^  = dU^  on  Sjj 


(10) 

(11) 


The  problem  thus  consists  in  determining  the  change  of  stress 
dd11  and  the  displacement  du^  in  such  a manner  that  the  boundary 
conditions  (10)  and  (11)  and  the  equation  of  equilibrium^)  are 
satisfied,  and  that  the  strain  change  competed  from  (8)  is  related 
to  the  given  values  of  0,  dO,  d^  and  h and  to  the  sought  stress 
change  dd^  by  means  of  Eqs.(4)  through  (7). 

The  first  extremum  principle  compares  the  actual  changes 
of  stress  and  strain,  dd^j  and  de^  , to  a fictitious  change  of 

4e  * 

stress  dd^j  and  the  corresponding  change  of  strain  de^j.  The 
stress  change  dd*j  is  supposed  to  satisfy  the  equation  of  equilib- 
rium and  the  boundary  condition  on  ST  , and  the  strain  change 
de*^  i?  associated  v/ith  ©,  dG,  d.^,  h,  and  dd!£,  by  means  of 
Eqs.(4)  through  (7),  but  need  not  be  derivable  from  a displace- 
ment field.  The  principle  of  virtual  work  then  furnishes  the 
equation 

j[(ddjj  « ddiWdeij  ] dV  = J[(dT*  - dTiJdUi  ] dSn  , (12) 

where  dT \ = dd  ^ j . 

The  integrand  of  the  left-hand  side  of  (12)  can  be 
transformed  as  follows; 
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♦ 

2(dd_».,  — ddjj)d£^^  — (ddj j dgjj  — dd.,^  de^^) 
XJ  XJ  x J •*■  J J — u 


— [ dd ^-j^dR-^j  ■•  de ^ ,| ) + de ^ ( dd ^ ^ — dd ^ ^ ) 1 • ( 13 ) 


13 


13 


On  the  righo-hand  side  of  this  equation,  introduce  the  com- 
ponents (4)  through  (?)  of  the.-  change  of  strain.  Since 

*Q  6 *q  c 

de  = de^  and  ds  ^ = de^j, 


♦ V-  jJcQ  jfcQ 

2(ddij  “ ddij)  deij  = ddi3(dei3  + dei3  + dei3^ 

" ddi.j(d£i.i  f de  ij  + dei3^  “ [dd*-j(de*|  + de*J  - de^  - dejj) 

+ ((36^+  dE^)(ddi;J  - dd *3  ) ] . (l4) 

The  bracket  in  (l4)  involves  only  elastic  and  plastic  changes 

* 

of  strain.  When  the  case  dd^j  = dd^j  is  excluded,  it  can  be 
shown  that  this  bracket  is  positive.  This  is  done  In  exactly 
the  same  manner  as  in  the  proof  of  the  extremum  princple  for 
isothermal  deformation  of  a work-hardening  plastic  material 
(see  Hill3,  pp. 63-64),  Thus,  the  left-hand  side  of  (14)  is 

4c 

smaller  than  the  first  two  terms  on  the  right-hand  side  unless  dd^ 
= dd^.  When  this  result  is  introduced  into  (12),  the  following 
relation  is  obtained* 

A /\ 

\ J [ddJjCdEij  + de *®  + dejj)  ] dV  - ] (dT^  dUi)dSu 


> 1 J CdcJij  (dejLj  + de^j  + dej.,)]  dV  - dUi)dSu  , (15) 

s|c 

where  the  equality  sign  holds  only  if  dd^  = dd^j.  The  relation 
(14)  establishes  a minimum  property  of  the  actual  change  of  state* 
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Before  the  second  extremum  principle  can  be  discussed, 
it  must  be  shown  that  Eqs.(4)  through  (7)  can  bo  transformed  so 
as  to  represent  dd^  as  function  of  de^.  We  note  first  that 
de^  and  de^  follow  immediately  from  the  data  of  the  considered 
boundary  value  problem.  Since  these  two  components  of  the  change 
of  strain  are  known,  giving  de^  is  equivalent  to  giving  de|.j 

P 0 C G 

+ de^j  = de^j  ” “ deij  • Tlie  Proof  that  the  sum  de^j 

+ de^  specifies  a unique  do^  then  proceeds  exactly  as  in  the 
case  where  de ^ ^ + de^j  represents  the  entire  change  of  strain 
(see  Hill3,  pp. 68-69). 

The  second  extremum  principle  compares  the  actual 
changes  of  strain  and  stress,  de^  and  dcs^  , to  a fictitious 


'iJ 


change  of  strain  dei;j  and  the  corresponding  change  of  stress 

★ * 

dd^  , The  strain  change  de^  is  supposed  to  be  derivable  from 

* 

a displacement  field  du^  that  satisfies  the  boundary  conditions 

]|c 

on  Sy  ; the  stress  change  dq^  is  associated  with  ©,  d©,  h, 


id 


and  de  . , by  means  of  F.q.(4)  through  (7),  but  need  not  satisfy 
the  equation  of  equilibrium  or  the  boundary  condition  on  S,p  . 

The  principle  of  virtural  work  then  furnishes  the  equation 

JtddijCde^  - de  ^ j ) ] dV  = [ [dTi(du*  - dut)  ] dSj.  (16) 

The  integrand  on  the  left-hand  side  of  (16)  equals 


/ *G  *P  e n . 

45  ij  (deij  + d£lj  - d£ij  - deiH 


(17) 
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and  does  not  involve  the  thermal  or  creep  effects;  it  can  there- 
fore he  transformed  in  exactly  the  same  manner  as  in  the  case 
where  the  sum  of  elastic  and  plastic  strain  changes  represents 
the  total  strain  change  (see  Hill^,  pp„ 65-66).  As  a result  of 
of  this  transformation,  the  following  relation  is  obtained: 

<f  (dTt  du*)dST  - J | [d<J*j(de*®  + de*jj)  ] dV 

P P 

< J (dTt  du.)dST  “ J d fddi^deij  + dEij^  dv»  (*8) 

* 

where  the  equality  sign  holds  only  if  de^  * deij*  The  relation 
(18)  establishes  a maximum  property  of  the  actual  change  of  state. 


